Nonlinearly-enhanced energy transport in many dimensional quantum chaos 



CO 

T— I 

o 

<N 



O 

o 



> 

00 

o 
in 

o 

m 



D. S. Brambila 1 ' 2 and A. Fratalocchifl 
1 PRIMALIGHT, Faculty of Electrical Engineering; Applied Mathematics and Computational Science, 
King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia 
2 Max- Born Institute, Max-Born-Stra(3e 2 A, 12489, Berlin, Germany 
(Dated: January 23, 2013) 

By employing a nonlinear quantum kicked rotor model, we investigate the transport of energy in 
multidimensional quantum chaos. Parallel numerical simulations and analytic theory demonstrate 
that the interplay between nonlinearity and Anderson localization establishes a perfectly classi- 
cal correspondence in the system, neglecting any quantum time reversal. The resulting dynamics 
exhibits a nonlinear ly-induced, enhanced transport of energy through soliton wave particles. 
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Anderson localization is a fundamental concept that, 
originally introduced in solid-state physics to describe 
conduction-insulator transitions in disordered crystals, 
has permeated several research areas and has become 
the subject of great research interest [H-[ll|. Theories 
and subsequent experiments demonstrated that disorder 
favors the formation of spatially localized states, which 
sustain diffusion breakdown and exponentially attenu- 
ated transmission in random media [1]. Although many 
properties of wave localization are now well understood, 
several fundamental questions remains. Perhaps one of 
the most intriguing problem is related to the transport of 
energy. Intuitively, one can expect that disorder — by fa- 
voring exponentially localized stated — arrests in general 
any propagation inside a noncrystalline medium. How- 
ever, the interplay between localization and disorder is 
nontrivial [H|, [l2| and under specific conditions random- 
ness can significantly enhance energy transport. In par- 
ticular, it has been observed that quasi-crystals with mul- 
tifractal eigenstates and/or material systems with tempo- 
ral fluctuations of the potential (or refractive index) , lead 
to anomalous diffusion in the phase space [IMsl. This 
originates counterintuitive dynamics including ultralow 
conductivities [13], as well as the formation of mobility 
edges even in one dimensional systems [H|- All these 
studies focused on specific geometries and linear mate- 
rials, while nothing is practically known about the role 
of nonlinearity in enhancing (or depleting) the transport 
of energy in disordered media. This problem acquires a 
strong fundamental character when refereed to the field 
of quantum localization. In this area, quantum-classical 
correspondences mediated by Anderson localization pos- 
sess many implications in the irreversible behavior of time 
reversible systems, which are at the basis of a long stand- 
ing physical dispute — i.e., the Loschmidt paradox [l7j — 
as well as many fascinating quantum phenomena such 
as the time reversal of classical irreversible systems and 
the quantum echo effect [Hi [H|- It has been argued, 
in particular, that microscopic chaos is at the basis of 
the irreversible entropy growth observed in classical sys- 
tems [20]. Time reversal, according to this interpreta- 



tion, is only possible at the quantum level TH, HI and 
sustained by Anderson localization, which breaks diffu- 
sive transport and suppresses the mixing ability of chaos 
[211 ]. However, when more dimensions are considered, 
numerical simulations predict that ergodicity is fully re- 
stored and diffusive transport settles is again, thus re- 
establishing the classical features of chaos and prevent- 
ing quantum time reversal [19]. Nevertheless, theoretical 
work reported to date considered only noninteracting sys- 
tems, characterized by linear equations of motion. The 
Loschmidt paradox, conversely, involved the use of inter- 
acting atoms, whose interplay in the mean field regime 
is accounted by short and/or long ranged nonlinear re- 
sponses 1 22- 24]. Besides that, as pointed out in the lit- 
erature [18| , atoms interactions are of crucial importance 
in quantum localization and diffusion. A key question 
therefore lies in understanding the role of nonlinearity in 
transporting energy in multidimensional quantum chaos. 
In this Letter, we theoretically investigate this problem 
by employing both numerical simulations and analytic 
techniques. To pursue a general theory, we here consider 
the following two dimensional model: 



• vV 



J drR(r f 



r)^(r , ) + /7^T(t) = 0, (1) 



with r = (x,y), V 2 = d 2 /dx 2 + d 2 /dy 2 , S T = En 6 ^ ~ 
nT) a periodic delta-function of period T, R a gen- 
eral nonlinear response and U(x,y) = 7(cosx + cosy) + 
ecos(x + y) a two dimensional periodic potential with 
strength defined by e and 7. Equation (pQ) defines a 
two dimensional, nonlinear quantum kicked rotor: for 
R = it reduces to the linear quantum kicked rotator 
[HI while for U = it corresponds to the 2D nonlinear 
Schrodinger equation (NLS), which represents a universal 
model of nonlinear waves in dispersive media [24J . In one 
dimension, conversely, Eq. (pQ) generalizes the nonlinear 
model investigated in [21] with classical chaos parame- 
ter K = 27T. Despite its deterministic nature, Eq. (pp) 
can be precisely mapped to the Anderson model with a 
random potential , and therefore furnishes a funda- 
mental model for studying energy transport in random 
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FIG. 1. (Color Online), (a)-(b) spatial density |^| 2 dis- 
tribution at (a) t = and at (b) t = 100T; (c) momentum 
diffusion (P) versus time in linear (dashed lines) and nonlin- 
ear (solid lines) conditions and for increasing coupling e. In 
the simulations we set a — 0.2, ojq — 0.3, A = 4 and K — 1.8. 



systems. The nonlinear response n = f drRir' — Y)ip(r r ) 
is modeled as a nonlocal term following a general diffu- 
sive nonlinearity (1 — <r 2 V 2 )n = |^| 2 , with nonlocality 
controlled by a. When a = 0, the system response is 
local with n = \ip\ 2 . For a ^ 0, conversely, the system 
nonlinearity becomes long ranged with kernel given by 
R(r) = ^ri^o(^), being Kq the modified Bessel function 
of second kind. Diffusive nonlinearities are particularly 
interesting in the context of nonlinear optics, as they can 
be easily accessed in liquids [26|, [27} , as well as in Bose- 
Einstein Condensates (BEC), where they generalize pre- 
viously investigated models [HI, [29| . 
We begin our theoretical analysis by calculating the mo- 
mentum diffusion (P) = (tp\^\*p) versus time, with p = 
V/i the momentum operator and (ip\f\ip) = J drf\ip\ 2 
the quantum average. Parallel numerical simulations are 
performed by a direct solution of (pQ) with an uncondi- 
tionally stable algorithm. In order for the field ip to ex- 
plore the periodic potential U, we here consider wave 
packets whose spatial extension Ar <C 2tt. Figure [T] sum- 
marizes our results obtained for a = 0.2, by launching 
at the input a gaussian beam ip = Ae~ x with waist 
co>o = 0-3 and amplitude A = 4 (Fig. [TJi). The stochas- 
tic parameter K has been set to K = 1.8 > if*, above 
the stochastization threshold K* & 0.97 where the lin- 
ear classical uncoupled rotor exhibits diffusive transport 
in momentum space [19|. For comparison, we also cal- 
culated the linear dynamics resulting from R = (Fig. 
[TJd dotted line). As seen from Fig. [TJd, the 2D nonlinear 
rotor behaves dramatically different with respect to its 
linear counterpart, demonstrating the strong role played 
by nonlinearity in the process. In particular, the linear 
system exhibits Anderson localization and diffusion sup- 
pression for e = (uncoupled condition), while for grow- 
ing e it shows a monotonically increasing sub-diffusion 
(Fig. [TJd). In the nonlinear regime, conversely, Anderson 
localization is suppressed even for e = 0, and the dynam- 
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FIG. 2. (Color Online). Positive Lyapunov exponent A 
versus coupling e, calculated for (a) Eqs. J4} and (b) Eq. J5}. 
In the simulation we set K — 5. 



ics shows an erratic, random-like behavior that does not 
manifest any simple monotonic increase for growing val- 
ues of e. These results are also significantly different from 
the nonlinear kicked rotor in one dimension [2l|, where 
nonlinearity was observed to induce localization effects. 
To theoretically understand this dynamics, we reduce the 
system to a nonlinear map modeling the nonlinear evo- 
lution of the ground state of Eq. (pQ). This analysis is 
justified by the observation that the spatial field profile, 
despite the chaotic motion, is not significantly altered in 
time (Fig. QJi,c). Due to the nonintegrability of the 2D 
NLS equation, we found analytical expressions by a vari- 
ational analysis [28|,[29|. In particular, we begin from the 
Lagrangian density C of Eq. (pp): 



i ( dip d>p* 



USrj 



M 2 

\ J dv'K(v-v'mv')f 



IWI 2 , (2) 



and study the ground state for U 

f^p e- r2/a2 



= by the following 
, defined by the power 



Gaussian ansatz: ip 

P = (^l^) an d waist a(t). By substituting the ansatz in 
([2]), after performing a variational derivative over a, we 
obtain a classical dynamics described by the following 
Hamiltonian Ti: 



V, 



2,71(7' 



: Z(a), (3) 



with Z(x) = e :r r(0, —x) and V acting as a potential for 
the one dimensional motion of a. The potential V has a 
bell shape profile that possesses a unique absolute mini- 
mum V(a*) for every combination of P and a. The fixed 
point a(0) = a* corresponds to a soliton wave of the 
system, which propagates in a translational fashion with 
fixed waist a(t) = a*, while different initial values lead to 
a breather [30] characterized by a periodic oscillation of a 
in time. When the kicks are turned on, for U ^ 0, the dy- 
namics of the ground state is perturbed by an addition of 
momentum p = (p x ,P y ), with a consequent translation of 
its center of mass. In order to model this dynamics, we 
considered the following general ansatz for the ground 
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FIG. 3. (Color Online), (a) Momentum diffusion (P) versus 
time calculated from Eq. (pQ) (solid lines), Eqs. (HJ)-© (dia- 
mond markers) and Eq. (pQ) in linear regime (dashed lines); 
(b) diffusion coefficient D versus coupling e. In the simula- 
tions we set a — 0.2 and K — 5. 



state evolution: = 



-(r-r ) / /a / +ip(r-r )/2T 



with 



p(t), a(t) and ro(t) = [xo(t),yo(t)] Lagrangian variables 
whose equations of motion, after an integration from nT 
to (n + 1)T, are found to be: 



Pn+l 

r n +i : 



= Pn - [7ng + e n usin(x + y ) 

r n + Pn+l 7 



(4) 



with 7 n = Ke"^, e n = 2eTe"^ L , g = [sinx , smy ], 
u = [1,1], f n = f( n T) and a n+ i calculated from the 
integration of the following equation: 

d^a _ dV 

dt 2 da T 

a 2 

x [7(cosxo + cos?/o) + 2e~~ e cos(xo + ?/o)]- (5) 

Equations (j4j) can be regarded as a variant of the four di- 
mensional standard map, which is randomized by time 
dependent coupling parameters j n and e n . The lat- 
ter depend on Eq. (j5j), which represents the motion 
of a two dimensional nonlinear kicked rotor. The sys- 
tem possesses an a dependent chaos parameter, given by 
K a = e- a2 / 8 7T. For K > K* , Eq. © is fully chaotic 
and can be regarded as an external noise source to Eqs. 
(|4|), increasing the mixing of the overall system [3l|. To 
highlight such a dynamics, we plot in Fig. [2^ and Fig. [2}d 
the positive Lyapunov exponent A [32] calculated for Eqs. 
(j4j) and Eq. ([5|), respectively. As seen in Fig. [2^, Eqs. (|4]) 
show a strong hyper chaotic behavior, with two positive 
Lyapunov exponents whose largest value grows linearly 
with e. Fig. Eb, conversely, displays the chaotic nature 
of wave packet extension a, whose Lyapunov coefficient A 
increases significantly fast (quadratically) with coupling. 
We investigate the diffusion in momentum p by observ- 
ing that above the stochastization threshold K > K*, 
the change in momentum Ap = p n +i — p n (xK becomes 
large compared to n. The classic position r n , which is 
taken modulo 27r, can be treated a random process, sta- 
tistically uncorrelated in time and uniformly distributed 
in [— 7r,7r]. The diffusion coefficient D is therefore evalu- 



ated as follows: 

D = ^^^) = ^!(sin^)( e -««/ 8 ) +2e 2 T 2 (e-°"/ 4 ) 

x (sm(x + y f) = ^(e" a »/ 8 ) + e 2 T 2 {e~ a ^ A ), (6) 

To evaluate the averages <e- a -/ 8 ) and {e~</ 4 ), we can 
consider a as a random variable (due to its chaotic motion 
in the phase-space), uniformly distributed between its 
oscillation extrema a min and a max : 



2 i 2 



erf 



V7TT 

~2Ar 



erf 



HO(4xA 3 ) (7) 



being A = a max — a min and having expanded the error 
functions up to second order, due to the smallness of 
their arguments a/r < 1. By substituting (O into (j6j), 
we obtain the diffusion coefficient, which reads as follows: 



D 



— + e 2 T 2 
4 



(8) 



Equation (j8j) allows to derive interesting properties for 
the nonlinear dynamics of Eq. (pQ). In particular, the 
quantum average (P) results from an hyperchaotic sys- 
tem described by a four dimensional standard map with 
random coefficients, and each realization manifests itself 
as a random walk in Fig. [TJd. The map diffusion rate 
is identical to the momentum diffusion of the classical 
linear rotor [l9| . hence, an additional average (in time 
or over an ensemble of input conditions) re-establishes a 
perfect classical correspondence for every coupling e > 0. 
It is worthwhile observing that the classical correspon- 
dence in the multidimensional linear quantum rotor is 
manifested only for very high coupling e, and in general 
the quantum diffusion (P) follows a fractional behavior 
with (P) oc (see e.g., [19] or Fig. [Ud dashed lines). 
As a result, the linear quantum rotor sub-diffuses at a 
slower rate than its classical counterpart. Conversely, 
Eq. (|Sj) predicts a perfect classical correspondence for 
every coupling e, which is re-established thanks to non- 
linear effects. In order to demonstrate this dynamics, we 
performed extensive numerical simulations from Eq. (pQ) 
and calculated the average diffusion through a quantum 
average followed by an average over different input con- 
ditions (P) = j D^{^W\^). Figure^ summarizes our 
results obtained for K = 5, a = 0.2 and by considering an 
initial wave packet composed by a Gaussian beam with 
waist ujo = 0.3 and amplitude A = 4. In complete agree- 
ment with Eqs. (|4|)-(|8|), we observe a diffusive behavior 
(P) cx t for every e > (Fig. [3^ solid lines and diamond 
markers), whose rate is significantly faster than the lin- 
ear subdiffusive dynamics (Fig. [3^ dashed lines). We 
can therefore conclude that nonlinearity favors the en- 
ergy transport in the system, increasing diffusion through 
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nonlinear wave-particles that are faster than their lin- 
ear counterparts. This result also highlights the intimate 
connection between the wave-particle aspects of nonlin- 
ear waves, whose quantum-classical characters cannot be 
singularly broken, but conversely emerge naturally after 
averaging over the corresponding degree of freedom. To 
further verify the scaling dependence predicted by Eq. 
dHJ), we calculated the diffusion coefficient D of the non- 
linear system for increasing e (Fig. [3b). In perfect agree- 
ment with our theory, we observe a quadratic behavior 
versus the coupling parameter e. 

In conclusion, motivated by the large interest in the 
study of energy transport in complex media, we inves- 
tigated the quantum-classical correspondences in many- 
dimensional quantum chaos. In particular, we em- 
ployed a two-dimensional, nonlinear quantum kicked ro- 
tor (NQKR) and study the role of nonlinearity in enhanc- 
ing or depleting energy diffusion and quantum-classical 
correspondences. We analytically tackled the problem by 
a variational analysis, reducing the dynamics to a four- 
dimensional standard map with random coefficients. In 
such an hyperchaotic system, a perfect classical corre- 
spondence is established by nonlinearity and an enhanced 
diffusion is observed due to solitons wave-particles, which 
are able to diffuse energy at a faster rate with respect 
to linear waves. These results show that quantum time 
reversal of classical irreversible systems is completely 
prevented in many dimensions, and demonstrate that 
nonlinearity can be effectively employed to increase the 
transport of energy in complex media. This work is ex- 
pected to stimulate further theory and experiments in 
the broad area dealing with quantum chaos and energy 
transport phenomena. 
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